Emergence of dissipative structures in current-carrying superconducting wires. 
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We discuss the emergence of a spontaneous temperature and critical current spatial modulation 
in current- carrying high temperature superconducting wire. The modulation of the critical current 
along the wire on a scale of 3 — 10 mm forces a fraction of the transport current to crisscross 
the resistive interface between the superconducting film and normal metal stabilizer attached to it. 
This generates additional heat that allows such a structure to be self sustainable. Stability and the 
conditions for experimental observation of this phenomenon are also discussed. 
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I. INTRODUCTION 

Propagation of electric current through a supercon- 
ducting substance is an inherently unstable process. If 
a small section of the superconductor turns normal, the 
heat generated by the current may lead to a thermal run- 
away process (quench). In practical (stabilized) super- 
conducting wires a normal conductor (stabilizer) is elec- 
trically and thermally coupled to the superconducting 
material. This allows the wire to recover from acciden- 
tal temperature perturbations or, if stabilization fails, to 
quench without irreparable damage [1]. Thus, a current- 
carrying stabilized superconducting wire is a bistable sys- 
tem with two spatially uniform stable modes of opera- 
tion: one is superconducting - low temperature and low 
(practically zero) dissipation, the other is normal - high 
temperature (greater than the critical temperature) and 
high dissipation. A transition from superconducting to 
normal mode initiated by a local temperature perturba- 
tion leads to a normal zone propagation (NZP) charac- 
terized by a certain speed. 

Here we will discuss the transition from superconduct- 
ing to normal mode of operation in the state-of-the art 
YBa2Cu2,07-:, {YBCO) coated conductors 0, (3. The 
main advantage of YBCO coated conductors over con- 
ventional low temperature superconducting wires - high 
operating temperature (65 — 77 K) - has an undesirable 
flip side: their heat capacity is high in comparison with 
that of low- Tc superconducting wires. As the result, the 
normal zone propagates slowly which significantly com- 
plicates protection of the devices, such as magnets or 
cables, from quench induced damage. We have investi- 
gated a possibility to increase the NZP speed by inserting 
a thin resistive layer (interface) between the supercon- 
ducting film and the stabilizer. There is always a very 
thin (a fraction of a micron) layer of material between the 
YBCO and copper that accounts for the resistance to the 
current exchange 0, [H . Conventional wisdom is that the 
interface resistance has to be as low as possible. If de- 
sired, as our findings suggest, it can be readily increased 
by various means. We have shown that the NZP speed 
can be substantially increased by increasing the interface 



resistance. This is the result of the greater amount of 
heat generated in the interface during current transfer 
from superconductor to the stabilizer [6]. Similar conclu- 
sions have been drawn in a large body of work, mostly 
theory, devoted to the effects of interfacial resistance on 
NZP in conventional superconductors [7, 8, 9]. 

The subject of this paper is a finding that when the 
interface resistance exceeds a certain threshold, the con- 
ventional scenario of transition between the supercon- 
ducting and normal modes of operation breaks down and 
a current-carrying superconducting wire may exhibit a 
much more complex behavior than it is commonly ex- 
pected. In addition to two uniform states - supercon- 
ducting and normal - there is also an anomalous nonuni- 
form, stable mode of operation characterized by a sponta- 
neously developing pattern of spatial temperature mod- 
ulation along the length of the wire with a peak temper- 
ature below the critical temperature Tc of the supercon- 
ductor. The temperature variation modulates the criti- 
cal current density, forming the static temperature and 
critical current density ripples (T-ripples or Jc-ripples). 
This is not a mesoscopic, but a macroscopic phenomenon, 
which can be classified as a dissipative structure - a result 
of the tendency for the spatially uniform physical systems 
driven away from thermal equilibrium to break the trans- 
lation symmetry and form steady macroscopic patterns, 
e.g. sand ripples[lO, Tl], Taylor- Couette flow, thermal 
convection, oscillatory chemical reactions, etc. [12]. The 
spatial scale of modulation is determined by the ther- 
mal diffusion length and in the case of coated conductors 
it is of the order of 0.3 — 1 cm, much greater than the 
magnetic field penetration or coherence lengths. 



II. PLANAR MODEL 

Coated conductors p|, [3| are manufactured in the form 
of a tape in which the superconducting YBCO film of 
about 1 /im thick is deposited on a buffered flexible metal 
substrate (e.g. Ni — W alloy, Hastelloy or stainless steel). 
A copper stabilizer is either soldered or electroplated on 
top of the YBCO film. Fig. 1. The typical width of such 



2 




transport the heat generated in the stabihzer. Thus, 



FIG. 1: A sketch of the cross-section of coated conductor (not 
to scale) . A thin superconducting film is sandwiched between 
the copper stabilizer (1) and metal substrate (2). The resistive 
interface between copper and YBCO and the insulating buffer 
between YBCO and substrate are not shown. 



a tape-like wire is 4 mm, the thickness, about evenly 
divided between the substrate and stabilizer, is close to 
100 fim and continuous pieces with uniform critical cur- 
rent density are manufactured in lengths over 100 m and 
up to 1 km.. 

In recent literature the normal zone propagation 
in coated conductors has been treated as a three- 
dimensional (3D) problem with the help of finite element 
analysis software, e.g. Ref.[l3]. In this section we show 
that the properties of the coated conductors allow the 
reduction of this nominally 3D problem to a 2D or ID 
problem without any significant loss of relevant physi- 
cal content. Computationally, such reduction makes the 
problem much more tractable, especially when it comes 
to the analysis of the effects of the interfacial resistance 
that introduces strong non-linearity into the partial dif- 
ferential equations describing NZP. 

The reduction of 3D equations of heat conduction in 
a thin tape-like composite wire to a 2D (planar) or ID 
(linear) model is possible if the variation of temperature 
across the thickness of the wire is negligible in comparison 
with its variation along the wire. The typical value of 
the current that coated conductor can carry is J ~ 200 — 
300 A /cm width. In the normal state the heat flux per 
unit area generated in the stabilizer 



(1) 



where pi and di are the resistivity and thickness of the 
copper stabilizer. At temperatures 77 — 90 K pi ^ 0.2 — 
0.3 pQ cm^ and for di = 40/im we get Jq ~ 1 W/cm'^. 
In order to remove this heat through the surface of the 
conductor a temperature gradient across the thickness of 
the stabilizer is required: 



KiAT 
di 



AT ' 



Ki 



(2) 



The thermal conductivity of copper Ki ~ 5W/cm 
so that the temperature variation across the thickness 
of stabilizer AT ~ 10~^ — 10~^ K. A similar estimate 
can be made for the substrate. Due to its relatively high 
resistance, the amount of Joule heat generated in the 
substrate itself is negligible [5]. The temperature gradi- 
ent across the thickness of the substrate is necessary to 
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do 
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diKi ' 
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The thermal conductivity of Hastelloy K2 ^ 7 x 
10~^W/cm K is substantially smaller than that of cop- 
per. Nevertheless, as long as the thermal flux jq is of the 
order of 1 — 10 W/cm^, the temperature variation across 
the substrate with thickness ^2 = 50 — 100 pm is still 
small, AT ~ 0.1 — 1 iC, in comparison with the variation 
of temperature along the conductor during NZP. This 
conclusion remains true for two other types of substrate, 
Ni — W alloy and stainless steel. 

The last structural element that needs to be exam- 
ined in terms of its effect on thermal conduction is the 
buffer between the YBCO film and the substrate. The 
buffer has a total thickness ^5 150 — 200 nm and con- 
sists of several layers of ceramic substances. If we take 
a representative value of the heat conductivity of Yttria 
stabilized Zirconia (YSZ) as ^ 1-5 x l^-'^W/cmK, 
the temperature drop across the buffer that is required 
to transfer the flux jq is 



AT, 



Kb ' 



(4) 



Again, as long as the thermal flux Jq is of the or- 
der of 1 — 10 W/cm^, the temperature drop AT5 ~ 
10~^ — 10~^i^. The thickness of the interface between 
copper stabilizer and YBCO is even smaller and it is also 
can be considered thermally transparent. Thus, the tem- 
perature of the coated conductors can be considered as 
dependent only on in-plane coordinates T(r), r = {x^y}. 

One of the specifics of coated conductors is that the 
superconducting film is much thinner than the stabilizer 
and the substrate. Therefore, it is natural to treat it as 
an infinitesimally thin layer, with negligible heat capac- 
ity capable, however, of carrying finite electric current 
density per unit width and generate a finite amount of 
heat per unit area. In a tape-like wire shown in Fig. 1 we 
place the superconducting film at z = 0, while the stabi- 
lizer occupies the volume < z < di, and the substrate 
is located at —d2 < z < 0. Since the estimates presented 
above have shown that in coated conductors the redis- 
tribution of thermal energy and equilibration of temper- 
ature in the z— direction is much faster than the heat 
propagation along the {x, y} plane, the two-dimensional 
equation of heat conduction can be obtained straightfor- 
wardly from energy conservation. The energy balance 
in a cylindrical volume of height di + d2 and arbitrary 
cross-section in the {x, y} plane leads to the following 
equation: 



dU 

~dt 



jQ,z\z=di JQ^ 



(5) 



Here U is the internal energy of the conductor per unit 
of its surface area, and 



dU_ _ dUdT _ dT 
'dt~df~dt^ 'di' 



(6) 
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where 



C = Cidi+C2d2 



(7) 



is the combined specific heat of the conductor. The con- 
tributions of the YBCO film, as well as that of buffer 
and the interface to the thermal mass of the conductor 
are negligible in comparison with that of the stabilizer 
and the substrate. 

The in-plane heat fiux is defined as 



JQ 



where 



K = Kidi^K2d2 



(8) 



(9) 



is the effective in-plane thermal conductivity of the con- 
ductor. The heat fiux from both, top and bottom sur- 
faces, in the right-hand side of Eq. (5), we will take in 
the Fourier form 



JQAz^-d^ =Ko{T-To). 



(10) 



Here Kq is the heat transfer coefficient across the insu- 
lation on the surface of the wire and Tq is the ambient 
temperature which, in the absence of losses, is the oper- 
ating temperature. In the situations when the conduc- 
tor is cooled with liquid coolant, the thermal ffux from 
the surface has substantially more complex temperature 
dependence[3,(8l]. But here we are interested in the effect 
of the interfacial resistance, so it make sense to consider 
all other contributing factors in their simplest form. The 
areal density of the internal heat sources in Eq. (5) 



Q 



q{z)dz 



(11) 



is the sum of all internal heat sources. Thus, the 2D 
(in-plane) heat conduction equations for the coated con- 
ductor takes form 



dT 



V ■ (KVT) =Q- 2Ko{T - To). 



(12) 



The redistribution of current between the supercon- 
ducting film and stabilizer is determined by the condition 
of charge conservation. 



V • J, + = 0; V -Ji-Jz 







(13) 



where Js and Ji are the linear density of current [A/cm] 
fiowing through the superconducting film and stabilizer 
respectively. The quench is a slow process in comparison 
to the time scale of charge redidtribution in a metal and 
the time derivative of the charge density in Eq. (13) can 
be neglected. The density of current fiowing across the 
interface between the stabilizer and superconductor 



J. 



Vi-Vs 
R 



(14) 



where Vi and Vs are the local electric potentials of 
the stabilizer and superconductor, respectively and 
R cw?] is the resistance of the unit area of the in- 
terface. 

There are three internal heat sources - originating in 
the stabilizer, interface, and superconductor respectively: 



1 

Pi 



El 



iVi - Vs. 
R 



-S{z)^Js-EsS{z). 



Here Ei = — VFi and Es = —Ws are the electric fields 
in the stabilizer and superconductor, respectively. We 
use the delta-functions to account for the fact that two 
of the heat sources are concentrated in the volume much 
thinner than either the stabilizer or substrate. The sub- 
strate contributes only to the thermal mass of the con- 
ductor due to its large resistance [5]. The integrated areal 
density of heat sources in Eq. (12) takes form 



Pi 



jVi - Vs? 
R 



Js - Eg. 



(15) 



Equations (12-15), supplemented by the nonlinear con- 
stituent relationship between current and electric field in 
the superconductor are sufficient to solve the 2D NZP 
problem. 



A. ID problem 

NZP can be treated as a one-dimensional problem 
when the width of the conductor is smaller or, at least, 
comparable to the thermal diffusion length, which we de- 
fine below. In this case the thermal equilibrium along the 
width of the conductor (in the direction) establishes 
more quickly than along the conductor, so that one can 
only consider the evolution of T(x, t), J(x, t) and Js(x, t) 
along the conductor. Since Ji -\- Js = J = const ^ the 
equations (13) reduce to one equation 



dJ^ ^ Vi-Vs 
dx 



R 



which can also be used in the form 



R 



= Es — El. 



(16) 



(17) 



The constituent relationship for a superconductor can be 
presented in many forms, the one most frequently used 
in literature is 



Es{Js) = Eq [ -J- 



(18) 



Here Jc(T) is the critical current and n is the expo- 
nent, usually large, n ~ 20 — 40. It is customary to 
take ^0 = 1/iV^/cm. For the stabilizer, the conventional 
Ohmic relationship will suffice: 



(19) 
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takes form 
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FIG. 2: The solid line is a sketch of the temperature de- 
pendent critical current density Jc. Temperature increases 
from left to right. The dashed line Js is the current density 
that flows through the superconductor. The dotted line is the 
constant total current density. The point Xc demarcates the 
boundary between the subcritical = and critical > 
sections. 



The one-dimensional version of Eqs. (12), (15), to- 
gether with Eq.(16) or (17), completely define the ID 
NZP problem. However, computationally, these two cou- 
pled equations still present a fairly formidable problem. 
Below we will introduce an approximation that allows us 
to reduce the problem to one equation that still preserves 
the relevant physics of the phenomenon. 



B. Reduction of ID NZP problem to one equation 

In good quality YBCO films the typical values of the 
exponent in Eq. (18) n ^ 20 - 40 [4]. We will de- 
fine the internal heat sources within the Bean model 
approximation [5, 9] which corresponds to the limit n 
oo. Correspondingly, as long as the current density Js 
is less than the temperature dependent critical current 
density Jc, the electric field and dissipation in the super- 
conductor is negligible. These we will call the subcritical 
sections of wire. In the critical sections of the conductor 
Js = Jc- In the normal state, T > Tc , no current flows 
through the YBCO film. 

Let us consider the situation shown in Fig. 2. The 
temperature rises and, correspondingly, the critical cur- 
rent density Jc{x) is falling from left to right. The point 
X = xi corresponds to the condition J = J^[xi). Well to 
the left of this point, in the subcritical section, all cur- 
rent flows through the superconductor. However, in the 
vicinity of xi the current starts to leak out of the su- 
perconductor into the stabilizer. Therefore, the critical 
section defined by the condition Js = Jc{xc) is located to 
the right of the point xi where J = Jc{xi). For x < Xc 
(the subcritical section) Es = and the condition of 
charge conservation, equivalent to Eq. (17) 



R 



dx'^ 



El — Es 



(20) 
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Pi J 
di 
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The solution of this equation 

Ji = Ae"(^-^^^ K = l/X= {pi/diRY^^ (22) 

has two unknowns, A and x^ These can be found by 
matching the in-plane current 



A = J -J, 



and out-of-plane current, see Eq.(13), 
nA 



dJc 
dx 



(23) 



(24) 



Let us introduce a dimensionless temperature and also 
assume for simplicity a linear dependence of the critical 
current on temperature |9|] 



T-Ti 
Tc-Ti 



; Jc{0) = J{l-0). 



(25) 



The current sharing temperature Ti is defined by the 
condition Jc(Ti) = J, so that T{xi) = Ti and, therefore, 
0{xi) = 0. Equations (23) and (24) take form 



f)f) 



(26) 



Thus, the current in the stabilizer in the subcritical sec- 
tion X < is given by 



Ji = JO{xc)e^^^-^^\ 
and in the critical section x > Xc 

Ji = J — Jc = JO. 



(27) 



(28) 



The location of the critical point Xc is determined by the 
condition 



ox 



(29) 



Here the ± sign corresponds, respectively, to either pos- 
itive or negative derivative dO/dx. The solution of the 
Eq. (29), 0{xc) = Oc, is always positive. Using the Taylor 
expansion 6{xc) ^ 0\xi){xc — xi) we get 



xi ~ A. 



(30) 



Hereafter a prime indicates a spatial derivative. 

Thus, in the subcritical region x < Xc the heat source, 
Eq. (15), is given by the sum of equal contributions from 
the stabilizer and interface 



Q = Pljl + R{J[f = j2(^(a;J)2e2-(--c)^ (31) 
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di 



di 
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In the critical sections of the wire where Jg = Jc ah 
three internal heat sources appear - originating in the sta- 
bilizer, interface and superconductor respectively. From 
Eq.(16) follows the relationship between the electric fields 
in the stabilizer and superconductor 



(32) 



The heat source in the critical sections x > Xc (Fig. 2) 
is given by 

Q = ^{J - Jcf + R(j:,f + (33) 

Jc (^(J - Jc) + i?Je') = ^ - Jc) + 

Expressed in terms of the dimensionless temperature, Eq. 
(25), the heat source in the critical section has the form 



di 



[e + \^{e')^-\\i-e)0"] 



(34) 



Finally, in the normal sections where T > Tc all current 
flows through the stabilizer and the heat source is given 

by 



Pi j2 



(35) 



In order to rewrite Eq.(12) in the standard dimension- 
less form and to reduce the large number of material con- 
stants to their relevant combinations, in addition to the 
deflnitions given by Eq.(25), we will measure the distance 
in units of thermal diffusion length It = (I^t/t)^^^, 
where Dt = K/C^ and time in units of 7~^, where the 
increment 7 = p\J^ j d\C 1ST determines the characteris- 
tic time required to warm an element of the conductor. 
Here AT = Tc- Ti 

Then, piece- wise deflned Eq. (12) at T > (l9 > 1), 
together with Eq. (35), takes form 



^-^ = l-.o{O-0o); e>l; 

to = 72 — ; T = 7*; ? = X It- 
Pi 



(36) 



In the critical sections with temperature Oc < < 1 
Eqs. (12) and (34) take the following form: 



dO d^O 



+ r{Oy - r(l - 0)0^' - Ko{0 - ^0); (37) 

er<o <1: 



Here Oc = 0{xc) is the floating boundary between the crit- 
ical and subcritical sections determined by the condition 
(29). The strength of the nonlinear terms is proportional 
to the interface resistance. 



A2 i? 

Lrp itQ 



Pill _ K{Te - Ti) 



di 



J2 



(38) 



Here A = {diR/ piY^'^ is the current transfer length, Eq. 
(22). 

Taking into account Eq. (31) we get for the subcritical 
regions 9 < On 



86 d^e 



2i9^e-2'^l^-^^l - A^o(6>-6>o); < 0^, (39) 



The equations (36)- (39), taken together, are piecewise 
deflned dimensionless version of the Eq.(12). Although it 
is more complex than the standard Kardar-Parisi-Zhang 
(KPZ) equation[14], they share a common feature - the 
competition between diffusion and the nonlinear growth 
term oc {0')'^. 

Even more important is that Eq. (37) may have nega- 
tive effective diffusion coefficient because it can be rewrit- 
ten as follows: 

1^ - (1 - r + rO)^ = e + r{e'f - ^,{9 - Oo); (40) 

(9c < 6> < 1; 

For r > 1 the effective nonlinear diffusion coefficient can 
be negative within the temperature interval Oc < < {r — 
l)/r. Therefore, the necessary condition for a negative 
diffusion coefficient is 



r > 



1 



>1, 



(41) 



where Oc is determined self-consistently by Eq. (29). 
Thus, only when the interfacial resistance, Eq. (38), ex- 
ceeds a certain threshold we can expect to see the pattern 
formation and other anomalous phenomena discussed be- 
low. 

To summarize, using the extreme limit of the con- 
stituent relationship (18), n ^ 00, we have reduced the 
system of the two equations, (12) and (13), to one piece- 
wise deflned equation (36-39) for temperature. The price 
we have paid for that is that the boundary between the 
critical {Es > 0) and subcritical {Es = 0) regions Oc is 
a floating one and has to be determined self consistently 
from Eq.(29). This presents a signiflcant challenge for the 
numerical solution. Our next step is to simplify the prob- 
lem by flxing the floating boundary at 6>c = 0. Comparing 
with the sketch in Fig. 2, this approximation means that 
for X < xi we take that all current flows through the 
superconductor, Jg = J, and for x > xi we take Js = Jc- 
This approximation becomes exact in the limit of small 
A <C /t, see Eq. (30), and we assume that even for X ^ It 
the qualitative results will provide a meaningful guidance 
to the physics of the phenomenon. 

The equation (36) which describes the normal section 
of the conductor does not change. In Eq. (39) we set 
Oc = 0. In Eq. (37) we need to ensure that all three 
heat sources deflned by Eq. (15) are positive. Specifl- 
cally it means that the electric fleld in the superconductor 
must have the same sign as the current. The direction 
of the current flow is taken to be positive. Therefore, 
Es = piJ/di(0 - X^O"), see Eq. (32), must be positive 
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FIG. 3: Sketch of the right-hand sides of Eqs. (36), (37), and 
(39) as a function of temperature in the Umit of uniform tem- 
perature. Sohd hues is the heat source, dashed hues indicate 
cooUng term for two values of coohng constant kq. 



or zero. Below we will use a step function to enforce this 
condition. Equation (37) will be used in the form 



r{e'f + (1 - e){e - rO")H{0 - rO") 

-^o{0-0o); 0<0<1; (42) 



Here the unit step function H{Es) ensures that the elec- 
tric field in the superconductor is either positive or zero. 
If the boundary between the critical and subcritical sec- 
tions were determined exactly, this condition would be 
met automatically. But since we set the boundary fixed 
at = 0, the step function is necessary to prevent the 
numerical solution to wander off into unphysical territory 
where the electric field Eg would have a direction oppo- 
site to current. In the range of temperature that meets 
the condition - rO" > Eq. (42) takes the form of Eq. 
(37) or (40). In the sections where H{Es) = the heat 
source is proportional to 6^ similar to that in Eq. (39). 



III. RESULTS 

The solutions of Eqs. (36), (39) and (42) with Oc = 
were obtained by two different methods in order to elim- 
inate the possibility of computational artifacts. We used 
a direct finite differences method and a high level soft- 
ware Mathematica. The temperature profiles ^(^, r) pre- 
sented below correspond to periodic boundary conditions 
and the initial condition in the form of a Gaussian in the 
center of the conductor 

0{^, 0)=eo + {\0o\ + a) exp{-eV2<^'}. (43) 

In all examples shown below we keep the width of the 
Gaussian on the order of the diffusion length, S = 
The dimensionless ambient (operating ) temperature is 
determined by the ratio of the transport current to the 



critical current at the operating temperature, Eq. (25), 

(44) 



^„ = i_:^<o. 



J 

The current sharing temperature is given by 



Ti = 



^0-1 



(45) 



The results shown below correspond to 6>o = — 1 (J = 
0.5Jc(To)). Correspondingly, the current sharing tem- 
perature T\ lies half-way between the operating and the 
critical temperature 



(46) 



The maximum temperature of the initial temperature 
profile, see Eqs. (25) and (43), is given by 



?;na^=Te + (a-l)(Te-Ti). 



(47) 



The results shown below correspond to the parameter 
a = 1.1. The respective maximum temperature of the 
initial temperature profile is slightly above Tc. The sta- 
bility of the conductor with respect to various values of 
the maximum temperature was reported in Ref.[6]. The 
system described by Eqs. (36), (39) and (42) has two sta- 
ble uniform {6' = 0" = 0) modes of operation, see Fig. 
3. One mode corresponds to the zero dissipation state 
of the conductor with temperature ^o- The other is a 
normal state with temperature 

^iv = ^o + ^o'>l; T^ = Ti + (^o + /^o')(^c-Ti) (48) 



The condition of bistability is 



1^0 < l^c 



1 



(49) 



For stronger cooling, /^o > ^c, Eq.(36) does not have 
a stable uniform solution and there is only one stable 
uniform mode of operation - the superconducting state 
with temperature ^o- 

The solutions of Eqs. (36), (39) and (42) with Oc = 
fall into several categories (scenarios) determined by 
the three dimensionless parameters {r, /^c^o}- Here we 
present a few scenarios that are important from both 
the basic physics and applications point of view. In 
all cases discussed below we take the transport cur- 
rent J = 0.5Jc(To), which corresponds to 6>o = —1 and 
tvc = 0-5. Figure 4 shows the 3D plots of the spatiotem- 
poral development of the initial Gaussian temperature 
profile for different values of the cooling constant tzo and 
the interface resistance above the threshold r = 1 (nega- 
tive effective diffusivity, Eq. (41)). 

In Fig. 4(a) = 0.1 (weak cooling. On = 9). This 
solution of Eqs. (36,39) and (42) describes a conven- 
tional bi-stable operation of the stabilized superconduct- 
ing wire. Additional heat generated in the interface 
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FIG. 4: Sequence of temperature profiles as they develop from the initial Gaussian perturbation visible in the foreground. All 
profiles correspond to the interface resistance r = 1.4, Eq. (38). The width of the initial temperature profile ^ = 1.4, Eq. (43). 
The range of the plots are < 7t < 50 and —40 < x//t < 40. (a) kq = 0.1; Conventional normal zone, but it propagates with 
greater speed than that at r ^ 1. (b) kq = 0.35; Conventional normal zone and temperature ripples emerge simultaneously, 
but propagate with different speed, (c) Cryostable condition {kq = 0.55 > k,c)- Only the temperature ripples are triggered. 



(KPZ-type growth) substantially increases the speed of 
normal zone propagation. The detailed analysis of the ef- 
fect of the increased interfacial resistance on propagation 
speed and wire stability is described in Ref.[6]. 

At stronger cooling, but still in the bistable regime 
(/^o = 0.35 < /^c), an anomalous third mode of conduc- 
tor operation emerges as shown in Fig. 4(b). The initial 
Gaussian profile gives rise to two distinct fronts propagat- 
ing with a different speed. A more rapidly propagating 
front leaves in its wake a metastable state characterized 
by the static spatial temperature modulation (T-ripples) 
with the peak temperature below Tc. The temperature 
modulation causes a modulation of the critical current 
density, which is why we may also call this pattern Jc- 
ripples. The Jc-ripples in turn are absorbed into the 
conventional normal state slowly propagating on top of 
them. 

When hzo > a^c, see Eq.(49) and Fig. 3, the cooling 
power at all temperatures is greater than the maximum 
power that can be dissipated in the stabilizer and, as the 
result, the uniform normal state is unstable. However, 
as Fig. 4(c) demonstrates, for r > 1 such cryostable, 
by conventional criterion, conductor is in fact a bi-stable 
system with the Jc-ripples being the second steady mode 
of operation. The initial temperature profile evolves into 
the Jc- ripples so that the average temperature of the 
conductor behind the front remains steady, close to Ti < 
Tc, instead of relaxing to the uniform superconducting 
state with temperature Tq. 



A. Properties of Jc-ripples 

In Fig. 5 the temperature profile of the Jc- ripples 
mode is shown. At a distance of a few diffusion lengths 
behind the propagating front the temperature modula- 
tion is stable and alternates in the range —0.8 < 6 < 0.26. 
Inspite of the overall symmetry of the problem, small 
computing errors and interpolation procedures lead to 
difference in the propagation speed of the two fronts and 



some irregularities in the profile. On the smaller scale, 
however, the profile is rather regular as shown in the in- 
set. The length of modulation is 2.7/t- 

A characteristic feature of these ripples is an apparent 
discontinuity of the first derivative dO/dx. This allows 
us to determine the value of the peak temperature Op = 
max{^(x)}. Integrating Eq. (40) over an infinitesimal 
interval that includes the temperature peak we get 



(1 - r + rOp) 



dO 



= 0. 



(50) 



peak 



Since the derivative is finite, this condition can be satis- 
fied by 



= (r- l)/r 



(51) 



This corresponds to zero value of the effective thermal 
diffusivity in Eq. (40). In numerical calculations the sec- 
ond derivative remains finite and the value of Op remains 
slightly below the limit given by (51). 

The heat generating regions {0 > 0) are spaced along 
the wire, segregated by the colder regions sustained by 
heat diffusion. One can compare the average power dis- 
sipation in the Jc-ripples mode and in the conventional 
mode. A section of the conductor between two minima 
of the temperature is thermally insulated from the rest 
of the system because the thermal flux vanishes at these 
points. Integrating Eq. (12) over the distance AL be- 
tween two troughs we get 

^jQdx = 2Ko{Tay - To); = ^ J T{x)dx 

(52) 

On the other hand, the power dissipation in the normal 
mode 

Qn = ^ = 2Ko{Tn - To) (53) 

Since T/v > Tav (see Eq. (48)), the average power dissi- 
pated by the Jc-ripples mode is smaller than that dissi- 
pated in the normal mode. In other words, in the regime 
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FIG. 5: Profile of the Jc-ripples (r = 1.4, kq = 0.4) at a 
certain moment of time, 7t = 25. The inset shows a magnified 
view of the same profile. The periodicity of the Jc-ripple 
is ~ 2.7 It- The peak temperature is Op ~ 0.24 — 0.26 in 
agreement with the exact value 1 — r~^,EQ. (51). 




FIG. 6: The initial condition gives rise to two propagating 
fronts that eventually merge. Here r = 1.5, a^o = 0.4, and = 
— 1. There is a domain boundary between the two merged 
fronts with slightly irregular temperature variation. 



where current is repeatedly exchanged between the su- 
perconductor and stabilizer the conductor behaves, at 
least in terms of losses, as a hyper conductor - a hypo- 
thetical substance with the resistance that is finite, but 
lower than that of the stabilizer, including the one made 
of copper. 

Figure (6) illustrates how two spreading domains of 
ripples, originating from different sources, interact. The 
initial condition is two Gaussians which do not overlap. 
Each of them gives rise to the spreading T-ripples. Once 
the spreading domains merge, a stable domain boundary 
characterized by a slightly irregular temperature varia- 
tion is formed. 




FIG. 7: Tri- stable operation which corresponds to = 0.4, 
r = 1.4, and Oq = —1. The initial profile evolves into Jc- 
ripples. A heater, Eq. (54), is switched on at jt = 20 and 
switched off at 'jt = 25. It triggers the transition to the 
normal state. The evolution of the system depicted in the 
Figure corresponds to 6 = 1.5, = 1.4 



B. External heat source and stability 

In order to study the stability of the ripples pattern 
we added a "heater" - an external heat source in the 
right-hand sides of Eqs. (36), (39) and (42) that can be 
switched on and off at the moments ti and t2 respectively: 

Qe,t = bexp{-e/2Sl}H{t - h)H{t2 - t). (54) 

Such a source simulates a typical experiment on normal 
zone propagation in which a small resistive heater is ther- 
mally anchored to the superconducting wire in order to 
trigger a transition to the normal state [l5]. 

Figure (7) illustrates the following scenario. Unlike in 
Fig. 4(b), the transition to the ripples mode is triggered 
by the initial condition without tripping the conductor 
into the normal state. Since the cooling constant = 
0.4 is below the cryostability criterion {hzo < = 0.5), 
the stable normal state with On = 1-5 is still possible. 
This indicates that there is a finite separation between 
the ranges of stability of the ripples mode and normal 
mode - one can be triggered without triggering another. 
At the moment = 20 the heater is turned on and 
turned off at 7^ = 25. A characteristic "rhino horn" of 
rising temperature triggers the secondary transition to 
the normal state. Thus, in this range of parameters, the 
conductor is a tri-stable, rather than a bi-stable system. 

An ability to trigger and observe the temperature rip- 
ples in an experiment depends, among other things, on 
stability and robustness of this pattern. A scenario shown 
in Fig. 7 has demonstrated a finite margin of stability of 
the ripples with respect to transition to the normal state. 
In Fig. 8 the effect of a " cold finger " - a negative power 
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FIG. 8: Effect of cold finger. Shown here is view "from 
above". The transition to Jc— ripples mode is triggered by 
the initial Gaussian profile. The negative power source col- 
lapses temperature in its vicinity, but does not destroy the 
ripple mode. Parameters here are the same as in Fig. 6. 

source is shown. Physically, this is equivalent to bring- 
ing the conductor in contact with a thermal mass with a 
temperature lower than Tq. The cold finger was modeled 
by a Gaussian, similar to that given by Eq. (54), but 
offset from the center and with a negative value of the 
parameter b. In Fig 8 b = —3, Si = 8, jti = 20, and 
7^2 = 35. After the application of the cold finger the 
temperature in its vicinity collapses well below the back- 
ground temperature = —1- The cold finger catches 
and holds the nearby propagating front, preventing it 
from expanding. Once the cold finger is lifted, the front 
resumes its propagation with the same speed as before. 
This demonstrates that the T-ripples mode is robust and 
recovers after severe disruptions. 



C. Stationary solutions 

The previous scenarios depicted in Figs. (4-8) show 
propagating fronts. Equations (36), (39) and (42) also 
have stationary solutions. A stable temperature profile 
can be established either from an initial perturbation or 
by the action of an external heater. Fig. 9(a) shows a 
simple bimodal structure similar to that shown in the in- 
set to Fig. (5). This bimodal soliton evolves from the 
initial Gaussian profile. A more complex five-peak struc- 
ture shown in Fig. 9(b) evolves under the action of the 
external heater, Eq. (54), with the same width as the 
initial condition in Fig. 9(a). In both cases the control 
parameter r = 1.2, Eq. (38), is the same, as are the 
cooling constant = 0.5 and the ambient temperature 
^0 = -1, Eq. (44). 

The stable solitons shown in Figs. 9(a,b) can be con- 




x/1t 

FIG. 9: View "from above", (a) A stationary bimodal so- 
lution that evolves from the initial Gaussian perturbation, 
(b) A more complex profile emerging under the action of the 
heater, Eq. (54). In both cases r = 1.2, a^o = 0.5, and 
Oo = -1. 



trasted with metastable - long lived perturbations shown 
in Fig. 10. In this scenarios the control parameter r 
is very close to the threshold r = 1. In Fig. 10(a) 
r = 1.0426. The perturbation evolving from the ini- 
tial Gaussian profile is quasistable and persists for about 
jt ~ 22 before dissipating. A very small increase of the 
interfacial resistance to r = 1.0427 yields a long lived per- 
turbation that eventually gives rise to normal zone prop- 
agation. This is consistent with our previous observa- 
tion that an increase of the interfacial resistance reduces 
cryost ability of the current carrying conductor: thermal 
perturbations that dissipate at lower values of resistance 
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FIG. 10: (a,b) Shown are long lived perturbations evolving 
from the initial Gaussian profile with ^ = 1.4. (a) The value 
of r = 1.0426 is very close to threshold r = 1 above which the 
effective diffusivity can be negative. The profile lingers for a 
long time, but eventually disappears, (b) At slightly greater 
r = 1.0427 the normal zone emerges and starts to propagate. 
In both cases kq = 0.35, = — 1 and ^ = 1.4. 

give rise to a NZP at higher resistance [g^. What is im- 
portant to realize comparing Figs. (9) and (10) is that 
it is the nonlinear KPZ term ~ (^')^ and the negative 
effective diffusivity in Eq. (37) that allow the isolated 
solitons like the ones in Fig. (9) to be stable and robust 
solutions over a finite range of parameters. 

D. Thermal oscillations 

Physical systems described by the parabolic partial dif- 
ferential equations do not support the propagations of 
waves. It is instructive, however, to investigate whether 
the nonlinearities may change that and allow a long range 
propagation of time-dependent perturbations originating 
from a localized source. The external power source, Eq. 
(54) was modified as follows 

Qe^t = 6exp{-^V2(^i'}(l + sin(2c^t))(l - e-^/"°). (55) 
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FIG. 11: (a) The ripple mode is triggered by the initial Gaus- 
sian profile. The pulsations of temperature with frequency 
2oo = 1.57 are caused by the heat source given by Eq.(55). 
Here r = 1.5, kq = 0.5, Oo = -1, S = Si = 1.4, 6=1. (b) 
The initial Gaussian profile gives rise to a bimodal station- 
ary soliton. The temperature pulsations are confined within 
the walls of the soliton. Here r = 1.2, a^o = 0.55, ^0 = — 1, 
6 = Si = 1.4, 6 = 0.3. 



This is equivalent to a resistive heater supplied with AC 
current with frequency u and switched on gradually over 
a period of time ~ tq. The results are presented in Fig. 
11. 

In Fig. 11(a) the ripple mode is triggered by an ini- 
tial temperature profile and the pulsating power source is 
gradually switched on with the time constant tq = 107~^. 
It is evident that the pulsations do not propagate outside 
the range of localization of the source. If there were a 
phenomenon of " thermal sound " it would manifest it- 
self as an undulation of the peak temperature lines that 
run along the time line. However, these lines are perfectly 
straight. 
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FIG. 12: The range of parameters {r, kq} specific for different 
regimes. The value of = —1 and the width of the initial 
perturbation ^ = 1.4 are the same for all combinations of 
{r,hio}. Area (1) corresponds to conventional NZP. In area 
(2) T— ripples and NZP coexist. In area (3) T— ripples are the 
second mode of operation. Area (4) corresponds to stationary 
solitons, and area (5) is the range of cryostability. 



Figure 11 (b) illustrates how the confinement of the 
thermal perturbations takes place. The control param- 
eter r = 1.2 and the cooling constant hzo = 0.55 are 
chosen to yield a bimodal stationary soliton similar to 
that in Fig. 9(a). The variable external power source 
creates temperature pulsations that remain strictly con- 
fined between the temperature peaks. This indicates that 
the singularities caused by the negative effective diffusiv- 
ity block the propagation of small perturbations through 
them. 



E. Phase diagram 

Figure 12 presents an approximate phase diagram that 
indicates the range of parameters {r, kq} where different 
regimes take place. This diagram corresponds to a par- 
ticular value of ^0 = — 1 which corresponds to transport 
current equal 50% of the critical current at ambient tem- 
perature To, see Eq. (44). It needs to be mentioned that 
the boundaries between the regimes depend not only on 
three parameters {r, hcq^Oq}^ but on the width of the ini- 
tial perturbation as well. This is obvious because the heat 
generated in the interface is determined by the derivative 
of temperature (the KPZ term in Eq. (42)). The dia- 
gram shown in Fig. 12 corresponds to the width of the 
Gaussian profile, Eq. (43), 5 = 1.4. 

In the range of parameters indicated as area #(1) in 
Fig. 12 the transition from superconducting mode to nor- 
mal takes the form shown in Fig. 4(a) - a conventional 
normal zone propagation, albeit with increased propaga- 
tion speed and reduced stability margins in comparison 
with that at low interfacial resistance. 

Area #(2) corresponds to simultaneous emergence of 



two propagating fronts as shown in Fig. 4(b). Area #(3) 
is a range of relatively strong cooling and interfacial re- 
sistance above the threshold r = 1. In this range of pa- 
rameters the conductor is bistable with T— ripples being 
the second stable mode of operation. Fig. 4(c). A strip 
#(4) corresponds to the stationary solitons. Fig. 9. Area 
#(5) is the range of cryostability - the temperature per- 
turbations dissipate without triggering a transition to the 
normal state. Figures 10 (a, b) illustrate the evolution of 
the initial temperature profile on both sides of the border 
between areas #(1) and #(5). 



IV. 



EXPERIMENTAL CONDITIONS 



In order to estimate the interface resistance required 
for the emergence of Jc- ripples we take the values of 
physical parameters from Ref.[4]. Let us take the op- 
erating temperature Tq = 67 K and Tc ~ 87 K. At 
J = (l/2)ji°\ Ti = (To + Tc)/2 = 77 K. The specific 
heat of copper and substrate at 77 K are, respectively, 
Ci ^ 1.7 J/cm^ K and C2 ~ 1.4 J/cm^ K. Let us take 
di — AO/j^m and ^2 = 50/im. Then, C = Cidi + C2<i2 ^ 

14.2 X 10~^J/ cw? K. The thermal conductivity of copper 
and substrate at 77 K are, respectively, Ki ^ bW/cm K 
and i^2 ~ 7 X l{)-^W/cm K. Thus, K = K^di + i^2^^2 ^ 

20.3 X lO-^W/K. The thermal diffusivity of the coated 
conductor Dt = K/C ^1.4 crn? / s. 

The self-field critical current density of the state of 
the art coated conductors at T = 67 i^T can be close to 
J^^ ~ 400 A/cm, and, correspondingly, the transport 
current density at 50% capacity can be J = 200 A/ cm. 
The sheet resistance of the stabilizer is pi/di ^ 0.5 x 
10~^Q. Taking AT = 10 we obtain the increment 7 = 
piJ'^/diCAT ^ 14 The thermal diffusion length 

It = (T^t/t)^^^ ~ 3 mm. The characteristic speed of 
NZP can be estimated as 



Ut = ItI ^ 4 cm/ s. 



(56) 



When the interfacial resistance is negligible, X <^ It^ the 
NZP speed is close to ut- This relatively slow propaga- 
tion speed complicates quench detection and quench pro- 
tection in coils made out of coated conductors. Increasing 
interfacial resistance leads to increasing NZP speed and 
reduced stability [6]. 

In order to observe the Jc- ripples the interface resis- 
tance has to exceed the characteristic value (Eq. 38)) 



Ro 



di 



5 X 10"^ Qcm^. 



(57) 



The interface resistance of the currently manufactured 
coated conductors is about 50 nQcm^[16]. Thus, two 
orders of magnitude increase in interface resistance is 
needed in order to create conditions under which the 
Jc- ripples may be observed. This can be readily ac- 
complished by modifying the stabilizer application pro- 
cedures, e.g. Ref.^l7j. 
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Taking into account the definition of the increment 7, 
the coohng constant a^q defined by Eq. (36) can be rewrit- 
ten as 



1^0 



2KoATdi _ 2Ko 



(58) 



To attain the value of hzo ^ 0.5, see Fig. 12, one needs to 
have the heat transfer coefficient 

Ko ^ ^ 5 X 10-2 W/cm^ K. (59) 

For a typical insulator like Kapton[4] the thermal con- 
ductivity at 77 K 



h^l.3x 10"^ W/cm K. 



(60) 



V. SUMMARY 



In summary, we have uncovered a novel pattern of su- 
perconducting wire operation in which the DC transport 
current is shared between the superconductor and stabi- 
lizer. This phenomenon manifests itself through sponta- 
neously developing temperature and critical current mod- 
ulation along the wire with the spatial scale of the order 
of 1 cm. Although the physics of this phenomenon is 
very different, the appearance of the temperature ripples 
is remarkably similar to that of the sand ripples. Most 
notably, this similarity extends to an apparent disconti- 
nuity of the slopes at the peaks [TO, [HI. 



Since 



do 



(61) 



where do is the thickness of insulation on the surface of 
the conductor, it is sufficient to have do ~ 200 — 300 fim 
to achieve the desired range of hco where one can expect 
to see the anomalous phenomena described above. 
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